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Using a novel, overdamped, oscillator flown aboard the Space Shuttle, we measured the 
viscosity of xenon near the liquid-vapor critical point in the frequency range 2 Hz </< 12 Hz. 
The measured viscosity divergence is characterized by the exponent z n = 0.0690 ± 0.0006, in 
agreement with the value z n = 0.067 ± 0.002 calculated from a two-loop perturbation expansion. 
Viscoelastic behavior was evident when t = (T - T c )/T c < 10~ 5 and dominant when t < 10“ 6 , 
further from T c than predicted. Viscoelastic behavior scales as /(/rwhere ris the fluctuation 
decay time. The measured value of A is 2.0 ± 0.3 times the result of a one-loop calculation. 
(Uncertainties stated are one standard uncertainty.) 
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As the liquid-vapor critical point is approached, the shear viscosity rj{C) measured in the 
limit of zero frequency diverges as £ n , where ^is the correlation length, which itself diverges 

on the critical isochore as r v . (Here, t = (T - T c )/T c , T c is the critical temperature, and v= 0.630.) 
Because the critical viscosity exponent z n = 0.069 is so small, it is very difficult to measure 
accurately on Earth in a pure fluid such as xenon. Close to the critical point, gravity causes 
stratification of the very compressible xenon and blunts the divergence of the viscosity in a 
manner that depends upon the height of the viscometer. Far from the critical point, the 
separation of the viscosity into critical and background contributions is uncertain; the critical 
contribution is a small fraction of the total and the separation depends sensitively upon the 
"crossover function," which is known only approximately. 

At non-zero frequencies f r/(£f) is complex in near-critical fluids, even at the sub-audio 
frequencies used here. Ordinarily, viscoelasticity (i.e. a partly elastic response to shear stress) is 
seen either at much higher frequencies or in complicated fluids such as polymer melts. 
Bhattacharjee and Ferrell [1] calculated the scaling function S(z) for near-critical viscoelasticity. 
Their result can be combined with the crossover function H(£) from Ref. [2] and the background 
viscosity tj 0 from Appendix C of Ref. [3] to obtain the prediction: 

nitf) = . ( 1 ) 

Here, the scaled frequency is defined by z = -inf; r, where the fluctuation-decay time is 
r = 6%Tj?/k B T c . We introduced the parameter A = 2.0 ± 0.3 into the argument of the scaling 
function S(z) to obtain agreement with our data. (All uncertainties stated are a single standard 
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uncertainty i.e. the coverage factor k = 1.) Accurate measurements of the scaling function 
require accurate measurements of ?]({;/) in the region where if is large while varying z by a large 
factor about unity. Prior measurements were unable to achieve these conditions [4,5]. 

We developed a novel viscometer, integrated it into the "Critical Viscosity of Xenon" 
(CVX) experiment package, and operated it aboard the Space Shuttle Mission STS-85. Figure 1 
shows that the Shuttle’s microgravity environment enabled CVX to make accurate viscosity 
measurements two decades in reduced temperature closer to T 0 than the best ground-based 
measurements. The data for the real part of the viscosity Re(^) lead to the result 
z^ = 0.0690 ± 0.0006, in agreement with the value z n = 0.067 ± 0.002 obtained from a recent two- 
loop perturbation expansion [6]. 

As CVX approached the critical temperature, viscoelasticity caused Re(rj) to depart from 
the asymptotic power-law divergence, first at high frequencies and then at low frequencies as 
shown in Fig. 2. In the same temperature range, Im(? 7 ) decreased from zero. (Im(^) is the ratio 
of the quadrature component of drag on the viscometer to the driving force.) This viscoelastic 
behavior was evident when t< 10‘ 5 and dominant at t < 10~ 6 , further from T c than predicted by 
Bhattacharjee and Ferrell’s one-loop calculation [1], The fitted value X = 2.0 ± 0.3 means that 
viscoelastic relaxation in xenon is 2 times slower than predicted. This factor is a challenge to the 
present one-loop theory; however, the challenge is comparable to one that Bhattacharjee and 
Ferrell met when scaling ultrasonic attenuation and dispersion in near-critical liquid mixtures. 
When they replaced the one-loop acoustic theory with a more sophisticated theory [7], the 
ultrasonic frequency scale changed by a factor of 1 .6. 

The unique CVX viscometer, shown schematically in Fig. 3, was a thin rectangular nickel 
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screen (19 mm long and 8 mm wide) composed of 0.03 mm wide electro-deposited wires spaced 
0.85 mm apart. A wire in the middle of the screen was extended and soldered to a supporting 
yoke, permitting the screen to oscillate in the xenon much like a child’s see-saw. These torsional 
oscillations were driven and detected by nearby electrodes. The torque applied to the screen was 
a repeated chirp with frequency components in the range 0.03 Hz </< 12.5 Hz. The 
corresponding viscous penetration length <5 = [//(/m/)]' 7 ' ranged from 800 pm to 40 pm. 

The screen had a complicated geometry; thus, the viscosity could not be calculated from 
its dimensions and its response to the applied torque. Instead, the drag on the screen was 
assumed to be proportional to f 2 B(<5) where the calibration function B(d) was determined from 
the viscosity data taken at ^cal = T c + 1 K while CVX was in orbit. The calibration procedure and 
tests of the assumed proportionality were described in our earlier work [8]. Because Im(^) = 0 
under the conditions of the calibration, we deduced Re( 7 ) and Im(^) from the magnitude and the 
phase of the screen’s response to the applied torques by analytic continuation of B(c 5) to complex 
arguments. 

We used data outside the viscoelastic region to check the calibration procedure in three 
ways. First, the values of Re(^) were examined at 2 Hz, 3 Hz, 5 Hz, 8 Hz, and 12 Hz for 
frequency dependence. None was found. Second, the values of Re(^) were compared with our 
previous data [3] taken on Earth with a high-(7, 7 mm high, oscillating-cylinder viscometer. 

They agreed throughout the range 3 * 10 4 < / < 5 * 1 O' 2 , where the data sets overlap and where 
the high-(7 data were not greatly affected by stratification. Third, the values of Im(^) were 
examined for departures from zero. The examination revealed departures from zero that were 
less than 0.002 Re(^) and were a weak function of log 10 (7’ - T cai ), suggesting that the screen was 
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influenced by an additional force not proportional to f 2 B{6). We applied a corresponding 
correction proportional to log, 0 (r-r cal ), and we extrapolated this correction to Im(^) inside the 
viscoelastic region. We checked the validity of this extrapolation by fitting A separately to Re(^) 
and Im(; 7 ) inside the viscoelastic region. The two values of A agreed within 4 %. (If we had 
assumed the correction to 1 01 ( 77 ) were proportional to (T -T c!ll ) raised to a small power such as z n , 
our results would be unchanged.) 

The maximum displacement of the screen was 0.03 mm and the maximum shear rate was 
9 s . The product (shear rate) x (fluctuation-decay time) was always sufficiently small that 
CVX did not encounter near-critical shear thinning [9]. The oscillating screen had a mass of 
only 1 mg; thus, its operation was relatively insensitive to the in-orbit vibrations of the Shuttle, 
and it survived the rigors of launch and landing with a large margin of safety. 

The CVX cell and thermostat were designed such that the xenon sample would be 
isothermal within 0.2 pK at T c ~ 290 K. Tests indicated that the thermostat’s performance was 
four times better than this [ 10 ], Calculations ignoring convection showed that the CVX sample 
was too large to attain thermal equilibrium very close to T c during the 1 1 -day flight. This led to 
the concern that once large density inhomogeneities were generated, for example by cooling the 
sample below T c , they would persist. To alleviate this concern, the sample’s temperature was 
programmed to minimize density inhomogeneities. In orbit, we discovered that large 
inhomogeneities decayed in approximately one hour at all temperatures above T c , as detected by 
the viscometer. Heat conduction along the metal components within the cell and electrophoretic 
convection of the sample contributed to this rapid decay. 

Key portions of the CVX data were taken as the temperature was ramped downward 
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through T c . The first ramp, at the rate - 1 pK/s, started at T c + 0.05 K. The second ramp was 20 
times slower; it started at T z + 0.003 K and had the rate -0.05 pK/s. Above T c , the data from the 
two ramps were indistinguishable; the viscosity increased monotonically, becoming complex and 
frequency dependent as T c was approached. Below T c , the data from the slower ramp indicated 
higher viscosities in the phase-separating sample. 

In Eq. (1), the background viscosity was described by the sum of analytic functions of 
the density pand the temperature: r/ 0 = t] w (T) + T) 0] (p). These functions were determined from 
the viscosity measured by others far from the critical point. (See Appendix C of Ref. [3].) For 
consistency with the measurements of Ref. [3] and the value p c = 1116 kg/m 3 , we set r} 0 (T c ,p c ) = 
(51.3 ± 0.4) pPa-s. The present results are not sensitive to this description of t) 0 [3,1 1]. 

The CVX viscosity data at 2 Hz, 3 Hz, 5 Hz, 8 Hz, and 12 Hz were used to determine 
five parameters. Two are the "universal" parameters z n and A; two are the wave vectors q c and 
q D that occur in the crossover function for xenon; and one is the value of T c on CVX’s 
temperature scale. When all of the data within the range 10 6 < t < 10 4 were fitted for these 
parameters, the results included: z n = 0.0690 ± 0.0006, q c £ 0 = 0.051 ± 0.007, and 
q D £ 0 = 0. 16 ± 0.05, where the standard uncertainties quoted allow for the correlations among the 
parameters. The same fit determined the frequency scale Az with a fractional standard 
uncertainty of 0.03. However, the fractional uncertainty of A itself is 0. 1 7 because of the larger 
uncertainty of the fluctuation decay-time amplitude r 0 given by 
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For xenon, r 0 = ( 1 . 1 6 ± 0. 1 7)* 1 O' 12 s and the uncertainty of r 0 is dominated by the uncertainly of 
xenon’s correlation length amplitude <f 0 = (0.184 ± 0.009) nm which we estimated [11] as the 
difference between the two measurements reported in [12]. 

The fit to the data more than 290 pK above T c determined T c with a standard uncertainty 
of 8 pK. It is reassuring that, for every frequency examined, both Re(^) and Im(^) showed 
extrema within 60 pK of the fitted value of T c . Equally reassuring, the fitted value of q c f 0 agrees 
with the value q c f 0 = 0.059 ± 0.004 determined from independently published data that do not 
include the near-critical viscosity [11], 

Fig. 2b shows the ratio Im( 7 )/Re( 7 ). This ratio is independent of the background 
viscosity tj 0 and the crossover function //(£); thus, these data are sensitive only to A, T c> and the 
scaling function. Fig. 2 demonstrates that the functional form S(Az ) describes both Re(^) and 
Im(77) in the ranges 10 6 < t < 10" 5 and 2 Hz < / <12 Hz, corresponding to 0.06 <Az< 33, 
provided that the factor A that sets the frequency scale has the value 2.0. 
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Captions for Figures. 


FIG. 1. Log-log plot of xenon's viscosity measured near the critical point. The asymptotic line 
has the slope y = z^v = 0.0435 deduced from the present microgravity data. Near T c , the CVX 
microgravity data [Re( 77) at 2 Hz] depart from the asymptotic line because of viscoelasticity. 

The two sets of ground data depart from the asymptotic line further from T c because the xenon 
stratifies in Earth’s gravity. 

FIG. 2. Xenon's viscosity at critical density measured at 2 Hz - 12 Hz. The solid curves resulted 
from fitting Eq. (1) to the data in the range 10" 6 < / < 10' 4 . (a) The real viscosity Re(^). Near 
t = 10 -5 , the data depart from the 0 Hz curve because of viscoelasticity, (b) The ratio 
Im(7)/Re(77). For clarity, the ratio data at frequencies above 2 Hz are displaced downward by 
integer multiples of 0.005; otherwise, they would coincide at / > 10' 5 . 

FIG. 3. Cutaway view of the CVX viscometer cell. The cylindrical volume occupied by the 
xenon was 38 mm long and 19 mm in diameter. The screen was oscillated by applying voltages 
to the four electrodes. 
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